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A recent polarized x-ray absorption experiment on the high temperature cuprate superconduc- 
tor P^S^CaC^Og+i indicates the presence of broken parity symmetry below the temperature, 
T*. where a pseudogap appears in photoemission. We critically analyze the x-ray data, and con- 
clude that a parity-breaking signal of the kind suggested is unlikely based on the crystal structures 
reported in the literature. Possible other origins of the observed dichroism signal are discussed. 
We propose x-ray scattering experiments that can be done in order to determine whether such 
alternative interpretations are valid or not. 

PACS numbers: 78.70.Ck, 78.70.Dm, 75.25.+Z, 74.72.Hs 



I. INTRODUCTION 

Twenty years since the discovery of high-temperature 
cuprate superconductivity, there is still no consensus on 
its origin. As the field has evolved, more and more at- 
tention has been directed to the pseudogap region of 
the phase diagram in undcrdoped compounds and the 
possible relation of this phase to the superconducting 
onei Time-reversal breaking has been predicted to occur 
in this pseudogap phase due to the presence of orbital 
currents^ and a subsequent experiment 3 - using angle- 
dependent dichroism in photoemission has claimed to de- 
tect this. However, this result has been challenged by 
others^ and an independent experimental verification of 
this would be highly desirable. 

Recently, Kubota et alfi performed Cu K edge circu- 
lar and linear dichroism x-ray experiments for under- 
doped Bi 2 Sr 2 CaCu208+ l; (Bi2212), claiming no time- 
reversal breaking of the kind predicted in Ref. Q exists, 
and that, on the contrary, a parity-breaking signal (but 
time-reversal even) is present with the same tempera- 
ture dependence as the photoemission dichroism signal, 
that was interpreted as x-ray natural circular dichroism 
(XNCD) as seen in other materials.— 

The aim of the present paper is to critically examine 
the conclusions of Kubota et al£ In particular, we find 
that the XNCD signal for the average^ Bi2212 crystal 
structure should be zero along all three crystallographic 
axes, therefore casting doubt on the original interpreta- 
tion of Ref. || To look into alternate explanations, we 
performed detailed numerical simulations aimed at ex- 
plaining the observed signal. At the basis of our study 
is the simple observation (see, e.g., Ref. H) that circu- 
lar dichroism in absorption can be generated either by 
a non-magnetic effect in the electric dipole-quadrupolc 
(E1-E2) channel (XNCD, a parity-breaking signal), or 
by a magnetic signal in the El-El channel (parity-even). 
Alternately, such a signal can be due to contamination 
from x-ray linear dichroism (XNLD). We propose x-ray 
experiments that could be used to investigate these mat- 



ters further. 

The structure of the paper is as follows. In Sec. II we 
show how symmetry constrains any possible XNCD sig- 
nal that would be observed in Bi2212. We also perform 
numerical simulations for XNCD assuming an alignment 
displaced from the c-axis, using several crystal structure 
refinements proposed in the literature. We also calculate 
the XNLD signal and comment whether XNLD contam- 
ination could be responsible for the observed signal. In 
Sec. Ill we calculate the x-ray magnetic circular dichro- 
ism (XMCD) signal at both the Cu K and L edges as- 
suming magnetic order on either the copper or oxygen 
sites. Finally, in Sec. IV, we draw some general conclu- 
sions from our work. 



II. NON-MAGNETIC CIRCULAR DICHROISM 
IN BI2212 

The structure of Bi2212 is strongly layered, with insu- 
lating BiO blocks intercalated between superconducting 
CuC>2 planes. Crystal structure refinements reveal the 
presence of an incommensurate modulation whose origin 
has been the subject of much debate. The presence of 
this modulation has complicated the determination of the 
average crystal structure. Two different average struc- 
tures have been proposed in the literature for Bi2212: 
Bbmb&i£ and Bb2b 11 ' 12 : 13 ' 14 , where b is the modulation 
direction for the superstructure. We follow the general 
convention in the cuprate literature and use a rotated 
basis compared to those in the International Tables for 
Crystallography 15 (respectively, Cccm, No. 66, and Ccc2, 
No. 37). In this way, the c-crystallographic direction is 
orthogonal to the CuC>2 planes. 

The Bbmb structure is globally centrosymmetric and, 
as such, does not admit a non-zero value for the parity- 
odd operator L ■ (e* x e)(f2 • fc), whose expectation value 
determines the XNCD signal (L, k, e and f2 are, re- 
spectively, the angular momentum operator, the direc- 
tion and polarization of the x-ray beam, and the toroidal 



2 



momentum operator, see, e.g., Refs. [§.16). Therefore, if 
the signal measured by Kubota et al^ were a true XNCD 
signal, this would imply a lower crystal symmetry than 
Bbmb. We note that most refinements in the literature 
suggest Bbmb, and this crystal structure is also con- 
sistent with recent photoemission dataJ^ that indicate 
the presence of both a glide plane and a mirror plane 
based on dipole selection rules. 

The other average structure that has been suggested 
by x-ray and neutron diffraction is Bb2b j 11 ' 12 i 13 ' 14 This 
space group is not centrosymmetric and, therefore, a par- 
ity breaking signal like that of XNCD is in principle al- 
lowed. However, not all wave vector directions are com- 
patible with the presence of a XNCD signal, as demon- 
strated below by symmetry considerations. In the last 
part of Section II, we numerically calculate the XNCD 
for a geometrical configuration allowing a signal - like 
A; 1 1 (101) - by means of the multiple-scattering subroutine 
in the FDMNES program^ 

In the context of this program, atomic potentials are 
generated using a local density approximation with a 
Hedin-Lundqvist form for the exchange-correlation en- 
ergy. These potentials are then used in a muffin tin 
approximation to calculate the resulting XANES sig- 
nal by considering multiple scattering of the photo- 
electron about the absorbing site within a one-electron 
approximation^ In the future, it would be desirable 
to repeat these calculations by using input from self- 
consistent band theory, as has recently been done for the 
Bbmb structure in regards to angle resolved photoemis- 
sion spectral 

In the Bb2b setting, Cu ions belong to sites of Wyck- 
off multiplicity 8d. These eight equivalent copper sites 
can be partitioned in two groups of four sites that are 
related by the vector (1/2,0,1/2). Within each group 
the four sites are related by the symmetry operations 
{E, Ci yi rh x , m z }, where E is the identity, C 2y is a two- 
fold axis around the b crystallographic axis, and rh x ( z ) 
is a mirror-symmetry plane orthogonal to the a(c)-axis. 
The absorption at the Cu K edge, expressed in Mbarn, 
can be calculated through the equations: 

- (±) =I>f\ (1) 

3=1 

of } = 4TT 2 ahu;J2 I <*i j) |0 (±) |*^> \ 2 S[hu; - (E n - E )] 

n 

(2) 

The operator (±) = e (±) •f(l + f£-r) in Eq. © is the 
usual matter-radiation interaction operator expanded up 
to dipole (El) and quadrupole (E2) terms, with the pho- 
ton polarization e and the wave vector k, where we label 
left- and right-handed polarization by the superscript ± . 

) i s the ground (excited) state of the crystal, 
and E (E n ) its energy. The sum in Eq. @ is extended 
over all the excited states of the system and Tllo is the 
energy of the incoming photon, with a the fine-structure 
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FIG. 1: (Color online) XANES signal for fc||(001) and XNCD 
signal for fc||(101) at the Cu K edge for three Bb2b crystal 
refinements, with a cluster radius of 4.9 A. The crystal struc- 
tures are Pet for Petricek et aZ.— , Glady for Gladyshevskii 
and Fliikigei^, and Kan for Kan and MossAS. The XNCD 
signals have been multiplied by the factors indicated. Each 
successive set of curves is displaced by 0.4 Mbarn. 

constant. Finally, the index j = 1,...,8 indicates the 
lattice site of the copper photoabsorbing atom in the 
unit cell. Eqs. |1| and ^ are the basis of the numer- 
ical calculations of the FDMNES program^ The eight 
contributions can be written as the sum of two equal 
parts coming from the two subsets of four ions related by 
the (1/2,0,1/2) translation. Within each subset, the four 
absorption contributions are related to one another by 
the symmetry operations 02 = C<i y a\, 03 = m x ai, and 
a 4 = rh z ai, implying that the total absorption is: 

a = 2{l+m z ){l + rh x )(Tx (3) 

The group from 05 to as is equivalent to the first group of 
four modulo a translation (this is the reason for the factor 
of two in Eq. ©). Notice that the symmetry operators in 
Eq. ([3]) are meant to operate just on the electronic part 
of the operator O in Eq. (J2J). 

In the case of circular dichroism, the signal is given by 
a = a + — a~ . If we suppose that no net magnetization is 
present in the material (we shall analyze the possibility 
of magnetism in Sec. Ill), then the dichroism is natural, 
i.e., necessarily coming from the interference E1-E2 con- 
tribution in Eq. ^ . In this case the signal is parity-odd, 
which implies that m z ( x ) = IC 2z (2x) — > -C 2z (2x) {I is 
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the inversion operator). Then Eq. ([3]) becomes 



(4) 



which implies that, of the possible five second-rank ten- 

(2) (2) 

sors involved in XNCD, only the term T| — T_{ sur- 
vives. To arrive at this result, we applied the usual op- 
erator rules on spherical tensors^! C2xTrn = and 

CzzTm = (-) m T, ( „ 2) . This in turn leads to a zero XNCD 
along the three crystallographic axes of the Bb2b crys- 
tal structure where, e.g., along the c-axis, the signal is 

(2) 

proportional to T . 

Therefore, even in the Bb2b crystal structure, the 
XNCD is exactly zero by symmetry when the wave vec- 
tor is directed along the c-axis (i.e., orthogonal to the 
Cu02 planes) as in the experiment of Ref. If! This has 
been further checked by numerical calculations with clus- 
ter radii up to 6.5 A, i.e., 93 atoms, centered on the Cu 
ion, based on the average crystal structures reported in 

Refs. |UHll!0. The only P° ssibilit y t0 J usti fy theo- 
retically the experimental evidence of circular dichroism 
of a non-magnetic nature is either by lowering of the or- 
thorhombic Bb2b symmetry, a misalignment of the k- 
direction with respect to the c-axis, or contamination 
from linear dichroism. We checked all of these possi- 
bilities. 

If we take into account the monoclinic supercell pro- 
posed in Ref. [H, with space group Cc, this implies a re- 
duction of the symmetry operations, with the loss of the 
two-fold screw axis. Nonetheless, the glide plane contain- 
ing the normal to the Cu02-planes is still present (rh x 
in Eq. [3]), which is responsible for the extinction rule of 
the quantity {^> n \L ■ (e* x e)(f2 ■ fc)|*„). Therefore the 
XNCD is again identically zero by symmetry, which we 
verified by direct numerical simulation of the supercell. 
Note that only a reduction to triclinic symmetry would 
allow for XNCD in the direction orthogonal to the Cu02 
planes. 

We also checked for the possibility of misalignment, as 
shown in Fig. [J by a direct calculation with fc||(101), 
corresponding to a tilting 9 ~ 10° with respect to the 
c-axis (note that the a lattice constant is ~ 5.4 A, and 
the c one ~ 30.9 A). We first remark that the calcu- 
lated XANES signal compares well to experiment (the 
XANES for | (001) and fc||(101) are identical on the 
scale of Fig. [T]). Despite this, the energy profile of the 
XNCD signal is different from the one reported by Kub- 
ota et alJi The main difference is the energy extension of 
the calculated dichroic signal, whose oscillations persist, 
though with decreased intensity, more than 50 eV above 
the edge itself. This characteristic is present for all four 
Bb2b refinements we have looked at (and as well for the 
monoclinic supercell, which we do not show) and is at 
variance with the experimental results of Ref. [f| where 
the dichroic signal is confined to a small energy range 
around the edge. 

Some comments on the calculations are in order. Each 
refinement gives rise to a different XNCD signal, and 
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FIG. 2: (Color online) XNCD signal, as in Fig. [TJ but for a 
CuCs cluster. 
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FIG. 3: (Color online) XNCD signal at the Cu K edge for 
fc||(101) as a function of the cluster radius (in A) for the crys- 
tal structure of Petricek et aZ.— . The signal has been mul- 
tiplied by 1000. Each successive curve is displaced by 0.06 
Mbarn. 



their intensities are quite different as well. They are 
found to be strongly dependent on the magnitude of the 
deviation of the atoms from their Bbmb positions, which 
differs significantly among the various Bb2b refinements. 
Moreover, the structure of Kan and Moss leads to a mani- 
festly different energy profile. This difference is seen even 
for CUO4 and CuO$ clusters (the results for a CuO^ clus- 
ter are shown in Fig. [2]) and has to do with the large de- 
parture of this particular structure from the Bbmb one. 
Although there has been some criticism in the literature 
concerning this particular refinement^ the point we wish 
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to make is that each refinement has a different XNCD 
signal, showing how sensitive this signal is to the actual 
crystal structure. We note that Fig. [1] was done for a 
cluster radius of 4.9 A, i.e., 37 atoms. In Fig. [3J we show 
results for the refinement of Petricek et al~ up to a ra- 
dius of 6.5 A (93 atoms), showing the development of ad- 
ditional structure in the energy profile as more and more 
atoms are included in the cluster. In the real system, the 
effective cluster radius is limited by the photoelectron 
escape depth, which is energy dependent^ 

We also remark that the magnitude of the XNCD sig- 
nal we calculated for a 10 degree misalignment is com- 
parable to that measured in Ref. 5- On the other hand, 
the signal goes as sin(2#), where 9 measures the displace- 
ment from the c-axis. We note that Kubota et al^ men- 
tion that their signal was insensitive to displacements 
from the c-axis of 5 degrees, which would argue against 
a misalignment given the strong angular dependence we 
predict. Moreover, we note that the size of the signal 
only depends on the projection of the k vector onto the 
a-c plane, i.e., a signal for fc||(lll) is equivalent to that 
for fc|| (101) . 

The above leads us to suspect that neither misalign- 
ment, nor symmetry reduction, are the basis of the signal 
detected in Ref. 0. We now turn to the third possibility 
for a non-magnetic signal, that due to intermixing of lin- 
ear dichroism. All x-ray beams at a synchrotron have a 
linear polarization component (Kubota et al£ mention 
the possibility of up to 5% of linear admixture). The 
resulting linear dichroism, which vanishes for a uniaxial 
crystal, can swamp the intrinsic XNCD signal for a biax- 
ial crystal where the a and b directions are inequivalent 
(like Bi2212). This was shown by Goulon et al. for a 
KTiOP04 crystal with the same point group symmetry 
(mm2) as Bi2212^ 

To analyze this further, we show the linear dichro- 
ism (the XANES signal for £J||(010) minus the one for 
E || (100)), calculated for the Bb2b refinement of Glady- 
shevskii et al.^* for various cluster radii, in Fig. 0J Simi- 
lar results have been obtained for the other crystal struc- 
tures, including the Bbmb refinement of Miles et aL— 

The energy profile, with a positive peak followed by a 
negative peak, and its location at the absorption edge, is 
very reminiscent of the data. Moreover, the size of the 
XNLD signal is large, meaning that only a few percent 
admixture is necessary to explain the size of the signal 
seen m Ref.i. One issue is that Kubota et al£ did report 
the existence of an XNLD signal, but also claim that it is 
temperature independent. This is somewhat puzzling, as 
there are significant changes of the lattice constants with 
temperature^ One obvious question would be why such 
an XNLD contamination would only appear below T*, 
though it should be remarked that there are anomalies 
in the superstructure periodicity near T*i2i A definitive 
test would be to rotate the sample under the beam, as 
any XNLD signal would vary as cos(2</>) where <f> is the in- 
plane angle relative to the &-axis. Any circular dichroism 
(XNCD or XMCD) is instead (j)- independent. 
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FIG. 4: (Color online) XNLD signal at the Cu K edge for 
fc 1 1 (001) as a function of the cluster radius (in A) for the crys- 
tal structure of Gladyshevskii et oZ.— . The XNLD signal has 
been multiplied by 50. The XANES curve has been displaced 
by 0.05 Mbarn. 



A final possibility would be a small energy shift be- 
tween the left and right circularly polarized beams. Dif- 
ferentiating the absorption edge in Fig. [T] would indeed 
lead to a signal similar to that seen in Ref. d (but with 
an enhanced positive peak relative to the negative peak). 
But such an energy shift is difficult to imagine with the 
particular experimental set up used. 

The observed dichroism signal as a function of energy 
is also reminiscent of that typically seen for magnetic 
circular dichroism: in this case the signal would be of 
El-El origin and its main features are indeed expected 
to be at the edge itself. In addition, the nature of the 
observed signal (a single sharp positive peak followed by 
a single sharp negative peak) is also much like a magnetic 
signal where the main features are expected to be more 
localized in energy starting from the rising edge of the 
absorption. Whether this possibility is realistic or not 
can only be determined by a quantitative calculation, 
which we offer in the next section. 



III. MAGNETIC DICHROISM IN BI2212 

Over the years, there have been several claims of pos- 
sible magnetic order in the pseudogap phase of cuprate 
superconductors. Recently, a magnetic signal at a (101) 
Bragg vector has been observed below T* for several un- 
derdoped YBa2Cu306+a; (YBCO) samples by polarized 
neutron diffraction^ The signal, corresponding to a mo- 
ment of order 0.05-0.1 /is, is not simple ferromagnetism 
as it was not observed at the (002) Bragg vector. Even 
more recently, a Kerr rotation below T* has been de- 
tected in underdoped YBCO corresponding to a net fer- 
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FIG. 5: (Color online) XMCD for fc || (001) at the Cu K edge 
with a ferromagnetic moment of 0.1 fiB along the c-axis at 
each copper site. The cluster radius is 4.1 A. The XMCD 
signal has been multiplied by 10000. The crystal structures 
are those of Fig. [1] 



romagnetic moment of 10~ 4 

This motivates us to consider the possibility of a mag- 
netic origin for the x-ray circular dichroism detected in 
Ref. H. We restate that in absorption (see, e.g., Ref. H), 
circular dichroism can be generated either by a non- 
magnetic effect in the E1-E2 channel (XNCD, a parity- 
breaking signal), or by a magnetic signal in the El-El 
channel (XMCD, parity-even). The first possibility, an- 
alyzed in the previous section, does not seem to be com- 
patible with the experimental results. In order to analyze 
the second possibility, we need to provide the lattice with 
a magnetic structure that has a net magnetization (oth- 
erwise, the XMCD is zero). In what follows, we shall sup- 
pose two magnetic distributions, the first with the mag- 
netic moments on the Cu sites, the second on the planar 
O sites. The numerical calculations are performed with 
the relativistic extension of the multiple-scattering pro- 
gram in the FDMNES code^ and provide results that 
are an extension of those of the previous section. 

The details of the calculations are as follows: we used 
again the average crystal structures discussed in Section 
II. For each magnetic configuration we employed clus- 
ters with radii ranging from 3.1 A (a CuOs-cluster), to 
4.9 A (37 atoms) around the Cu photoabsorbing ion. In 
the first set of calculations, shown in Figs. [5] and [HJ we 
built the input potential from a magnetic configuration of 
4.55 3df electrons and 4.45 3<i| electrons (i.e., a moment 
of 0.1 fiB per copper site). In the second set of calcula- 
tions, shown in Fig. we built the input potential from 
a magnetic configuration of 2.05 2p^ electrons and 1.95 
2pi electrons (i.e., a moment of 0.1 /is per planar oxygen 
site). The following results are noteworthy: 

a) Differently from the XNCD calculations shown in 
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FIG. 6: (Color online) XMCD as in Fig. [5] but with a cluster 
radius of 3.1 A, i.e., a CuOs cluster. 
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FIG. 7: (Color online) XMCD for fe||(001) at the Cu K edge 
with a ferromagnetic moment of 0.1 fiB at each planar oxygen 
site. The XMCD signal has been multiplied by 5000. The 
XANES signal has been displaced by 0.02 Mbarn. 



Fig. [I] all crystal structures give basically the same 
XMCD spectra. The reason for this behavior may be 
related to the fact that XMCD, when x-rays are orthog- 
onal to the Cu02 planes, mainly depends on the in-plane 
magnetisation density and on the in-plane crystal struc- 
ture, which is quite similar for the various refinements. 

b) There is a more marked dependence on the cluster 
radius compared to the XNCD, as shown by the compar- 
ison of Fig. [5] and Fig. [6j The calculations with a radius 
bigger than 4.1 A (6 Cu, 9 O, 4 Sr, and 4 Ca), above the 
pre-edge energy, are basically equivalent to those shown 
in Fig. [5J with a positive peak at the edge energy, fol- 
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FIG. 8: (Color online) XMCD for fe||(001) at the Cu L 2 , 3 
edges with a ferromagnetic moment of 0.1 fiB along the c- 
axis at each copper site. The cluster radius is 4.1 A. The 
XMCD signal has been multiplied by 10. 



lowed by a double negative peak, the latter at variance 
with the experimental results. On the contrary, the en- 
ergy shape obtained for a radius of 3.1 A is very close to 
the experimental one, with a single negative peak after 
the sharp positive one, with relatively good agreement 
in the energy position and width. We could be tempted 
to suppose, therefore, that the virtual photoelectron has 
a very small mean free path before decaying and it is 
sensitive just to the nearest neighbor oxygens. Indeed 
we checked that an identical XMCD profile is obtained 
with just the in-plane Cu04-cluster. On the other hand, 
the size of the signal we calculate is about an order of 
magnitude smaller than that seen in Ref. [H. Since the 
XMCD signal is proportional to the moment, then we 
would need a moment of ~1 Bohr magneton per copper 
to have a comparable signal. Such a huge moment would 
have been observed previously by neutron scattering if 
it existed. Of course we cannot exclude that spurious 
effects, such as strain fields, could have influenced the 
measurement. 

c) The energy profile in the case of magnetization at 
the oxygen sites is not much different from the cop- 
per case, except for the deeper negative peak around E 
~ 8.994 keV, as shown in Fig. [7] Also in this case the dif- 
ferent crystal structure refinements give basically equiv- 
alent results, as again the Cu02 planes are practically 
equivalent in the various cases. Note that the relative 
intensity is equivalent to the copper case, as 0.1 Bohr 
magnetons per planar oxygen corresponds to 0.2 Bohr 
magnetons per Cu02 cell (note we multiply by 5000 in 
Fig. [7] as compared to 10000 in Fig.©. 

We also performed simulations for the Cu 1^3 edges 2 ^ 
for the magnetic configuration corresponding to Fig. O 



as shown in Fig. [8l which can be compared with future 
experimental investigations in order to confirm whether 
or not a net magnetization exists in this compound. 

Finally, we remark that the dependence of the XMCD 
signal on the tilting angle 9 (i.e., the displacement of 
the photon wave vector from the c-axis) goes like cos(0) 
and therefore the signal is not very sensitive to small 
displacements of 5 degrees, as noted by Kubota et al& 
This different angular dependence from the XNCD signal 
suggests a relatively easy way to unravel the question ex- 
perimentally: it is sufficient to measure the 9 (azimuthal) 
dependence of the signal, noting that any XNLD contam- 
ination would be tested by the <f> (polar) dependence of 
the signal. 



IV. CONCLUSIONS 

In our opinion, the experiments of Kubota et al& have 
raised more questions than they have answered. Al- 
though not treated in our paper, we believe that their 
results from the measurement of non-reciprocal linear 
dichroism are at this stage not conclusive, as only one di- 
rection for the toroidal moment has been investigated, of 
the two possible suggested by the orbital current pattern 
of Varma. 2 The analysis performed in Section II showed, 
moreover, that the claimed XNCD signal is probably un- 
justified. In fact, even though the space group Bb2b 
is non-centrosymmetric, XNCD is absent by symmetry 
when the x-ray wave vector is chosen orthogonal to the 
Cu02 planeSj as in the experimental measurement geom- 
etry of Ref. [5J. The same extinction rule survives for the 
monoclinic supercell structure refined in Ref. [l3l More- 
over, in both cases, it seems hard to mantain the hy- 
pothesis of misalignment, due to the experimental local- 
ization of the energy profile around the main absorption 
edge, which is absent in the calculations. We also note 
the difference of XNCD from the photoemission dichro- 
ism results of Kaminski et al£ A direct comparison is 
however not immediate, as the former represents a q- 
integrated version of the latter (here q is the solid angle in 
the space of the photoelectron wave- vector, see, e.g., Ref. 
8). A more likely explanation is an XNLD contamination 
(Fig.dJ), but then the challenge is to understand why such 
an effect would only exist below T*. We note that an op- 
tics experiment for an optimal doped Bi2212 sample has 
seen a change in linear birefrigence below T c , which was 
accompanied by a non-zero circular birefringence^ In 
addition, both Bi2223 29 and the Fe analogue of Bi2212^ 
exhibit supercells with 222 space groups which would al- 
low for dichroism. So, it is conceivable that there is a 
subtle structural transition associated with T* which we 
suggest could be looked for by diffraction experiments. 

A final comment about the physical quantities detected 
by x-ray circular dichroism, either natural or magnetic, 
is in order. At the K edge of transition metal oxides, 
XMCD in the El-El channel, at the excitation energy 
£ = Tilo — -Eedge, gives information on the expectation 
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value of L • (e* x e) for the excited states at the energy £ . 
The orbital angular momentum is either induced from a 
spin moment via spin-orbit coupling (as calculated here) 
or directly by an orbital current (as in the scenario advo- 
cated in Ref . • No direct spin information is available 
at the K edge and therefore an XMCD measurement is 
not directly related to the ground state magnetic moment 
as for the L edge. Moreover, the observed states are 
those with p-like angular momentum projection on the 
photoabsorbing Cu ion, which are extended and there- 
fore mainly sensitive to the influence of the oxygen atoms 
surrounding the Cu site. In this case, the main contribu- 
tion to the XMCD energy profile is expected in an energy 
range of 10-20 eV from the main edge to the first shoul- 
der in the XANES spectrum, as found in Ref. and in 
our own XMCD calculations. 

Although the results of Section III are in principle con- 
sistent with Ref. 0, the size of the ferromagnetic moment 
necessary to get a signal of the magnitude seen by exper- 
iment, ~1 Bohr magneton, is excessive. If such a large 
ferromagnetic moment existed, it would have surely been 
seen by neutron scattering. From this point of view, ex- 
periments performed at the Cu L edge and O K edge 
would be desirable as they are more sensitive to the pres- 
ence of a magnetic moment. 

Finally, we would like to remark that XNCD along 
the c-axis is insensitive to orbital currents. These latter, 
confined to the Cu02 planes, develop a parity breaking 
characterized by a toroidal moment (f2) within the Cu02 



planes. The XNCD experiment of Ref. [H would only be 
sensitive to the projection of the toroidal moment out of 
this plane (i.e., along the direction of the x-ray wavevec- 
tor). Therefore, if performed as stated, it cannot tell us 
about any possible orbital current order. 

To conclude, we believe that the various interpreta- 
tions, XNCD, XNLD, or XMCD, have their drawbacks, 
and therefore the origin of the experimental signal of 
Ref. Ifa is still open. In this sense, further experimen- 
tal checks of the energy extension of the dichroic signal 
would be highly desirable. Based on our results, the most 
stringent experimental test on the physical origin of the 
signal would come from the measurement of the depen- 
dence on the tilting (azimuthal) angle 6, due to the dif- 
ferent dependences of XNCD and XMCD, as well as the 
dependence on the in-plane (polar) angle 0, which would 
test for any possible XNLD contamination. 
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